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1 Introduction: connection of fluctuations and nonextensivity 

Both the notion of fluctuations and that of nonextensivity are nowdays widely 
known, albeit mostly in the fields of research only indirectly connected with 
particle physics. Both turns out to be very fruitful and interesting and this 
is cleaxbf demonstrated by all other lectures given at this workshop (see espe- 
cially ffl). 

This lecture will be devoted to the case in which evident nonextensivity 
of some expressions originate in intrinsic fluctuations in the system under 
consideration (the origin of which is usually not yet fully understood)^. The 
best introduction to this problem is provided by the observation that in some 
cosmic ray data-Hike depth distribution of starting points of cascades in Pamir 
lead chamber E3) one observes clear deviations from the naively expected 
exponential distributions of some variables which are evidently better fitted 
by the power-like formulas: 
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Her£.iV denotes the number of counts at depth h (cf. El for details). Whereas 



in cij we have proposed as explanation a possible fluctuations of the mean 

free path A in cq. (|l|) characterised by relative variance ui — - — — - > 

0.2, in i the same data were fitted by power-like (Levy type) formula as 
above keeping A fixed and setting q = 1.3. In this way we have learned 



a Our encounter with this problem is presented in works 
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about Tsallis statistics and Tsallis nonextensive entropy and distributions^. 
By closer inspection of the above example we have been able to propose a new 
physical meaning of the nonextensivity parameter q, as a measure of intrinsic 
fluctuations in the system □ Fluctuations are_.therefore proposed as a new 
source of nonextensivity which should be added £□ to the previously known and 
listed in literature sources of the nonextensivity (like long-range correlations, 
memory effects or fractal-like structure of the corresponding phases space □) . 

To demonstrate this conjecture let us notice that for q > 1 case, where 
e G (0, oo), one can write a kind of Mellin transform (here a = Q: 




where f q> x(l/X) is given by the following gamma distribution: 




with /i = aAo and with mean value and variation in the form: 




where / 9 <i(l/A) is given by the same gamma distribution as in but this 
time with a — ► a' and /i = fj,(e) — a'Xo — e. Contrary to the q > 1 case, this 
time the fluctuations depend on the value of the variable in question, i.e., the 
mean value and variance are now both e-dependent: 



'See Tsallis lecture i and references therein (cf. also B0) for detailed necessary information 
concerning Tsallis statistics and the non-extensivity. 



crete 'q: submitted to World Scientific on February 1, 2008 2 



However, in both cases the relative variances 
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remain e-independent and depend only on the parameter q |j. It means there- 
fore that oO (at least for the fluctuations distributed according to gamma 
distribution) 

L=exp (~£) ^ =exp «(-^)=( exp (-f)>' (s) 

with q = 1 ± uj for q > 1 (+) and q < 1 (— ), i.e., there is connection between 
the measure of fluctuations lu and the measure of nonextensivity q (it has been 
confirmed recently in EJ) . 

2 Where are the fluctuations coming from? 

2. 1 Generalities 

The proposed interpretation of the parameter q leads immediately to the fol- 
lowing question: why and under what circumstances it is the gamma distribu- 
tion that describes fluctuations of the parameter X ? To address this qjj^tion 
let us write the usual Langevin equation for the stochastic variable A 



dX 
~dt 



- + m 

T 



A = (f> = const > 0. (9) 



with damping constant t and with source term <p, different for the two cases 
considered, namely: 

, , 1 /\ £ \ , , , A 

<P = <Pq<\ = 



X j j whereas 4> — 0g>i = — • (10) 



For the usual stochastic processes defined by the white gaussian noise form 
of £ (t) one obtains the following Fokker-Plank equation for the distribution 



c Notice that, with increasing a or al(i.e., for q — t 1) both variances (Jij) decrease and 
asymptotically gamma distributions becomes a delta function, /(1/A) = <5(A — Aq). 
d \t means that ensemble mean (f(t)) = and correlator (for sufficiently fast changes) 
+ At)) = 2 D (5(Ai). Constants r and D define, respectively, the mean time for 
changes and their variance by means of the following conditions: (A(t)) = Aq exp (— |r) 
and (A 2 (t = oo)) = |Dt, Thermodynamical equilibrium is assumed here (i.e., t >> r, 
in which case the influence of the initial condition vanishes and the mean squared of A has 
value corresponding to the state of equilibrium). 
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function of the variable A 

*a + ( u) 

where the intensity coefficients K\ t 2 are defined by eq.(^) and are equal toSH: 
KAX) = <j) - - + D\ and K 2 (X) = 2DX 2 . (12) 

T 

From it we get the following expression for the distribution function of the 
variable A: 



/(A) = Jfc(A) 6XP 
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which is, indeed, a gamma distribution ([|) in variable 1/A, with the constant c 
defined by the normalization condition, J Q °° d(l/A)/(l/A) = 1, and depending 
on two parameters: fi(s) = ^sdEi an d a q — with <j) q = 4> q ^x,q<x an d 
a q = (a, a') for, respectively, q > 1 and q < 1. This means that we have 
obtained eqs.(^) with ui — and, therefore, the parameter of nonextensivity 
q is given by the parameter D and by the damping constant r describing the 
white noise. 

2.2 Temperature fluctuations 

The above discussion rests on the stochastic equation (||) . Therefore the pre- 
viously asked question is not yet fully answered but can be refrazed in the 
following way: can one point the possible physical situation where such fluc- 
tuations could be present in the realm of particle physics ? Our proposition to 
answer it is to identify A = T, i.e., to concentrate on the possibility of fluctu- 
Mions of temperature widely discussed in the literature in different contexts 
t3. In all cases of interest to us temperature T is variable encountered in sta- 
tistical descriptions of multiparticle collision processes t£l. Our reasoning goes 
then as follows: Suppose that we have a thermodynamic system, in a small 
(mentally separated) part of which the temperature fluctuates with AT ~ T. 
Let \{t) describe stochastic changes of the temperature in time. If the mean 
temperature of the system is (T) = Tq then, as result of fluctuations in some 
small selected region, the actual temperature equals T 1 = Tq — r£(t)T and 
the inevitable exchange of heat between this selected region and the rest of 
the system leads to the process of equilibration of the temperature which is 
described by the following equation of the type of eq. (^|) E3 : 

^-\{T> -T) + Q q = (14) 
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(here 51 g< i = and f2 ? >i = 0). In this way we have demonstrated that eq. 
(^) can, indeed, be used in statistical models of multiparticle production as 
far as they allow the g-statistic as their basis. 

The q < 1 case requires a few additional words of explanation. Notice 
the presence of the internal heat source, which is dissipating (or transfering) 
energy from the region where (due to fluctuations) the temperature T is higher 
to colder ones. It could be any kind of convection-type flow of energy; for 
example, it could be connected with emission of particles from that region. 
The heat release given by e/(ra') depends on e (but it is only a part of e 
that is released) . In the case of such energy release (connected with emission 
of particles) there is additional cooling of the whole system. If this process 
is sufficiently fast, it could happen that there is no way to reach a stationary 
distribution of temperature (because the transfer of heat from the outside can 
be not sufficient for the development of the state of equilibrium). On the 
other hand (albeit this is not our case here) in the reverse process we could 
face the " heat explosion" situation (which could happen if the velocity of the 
exothermic burning reaction grows sufficiently fast; in this case because of 
nonexistence of stationary distribution we have fast nonstationary heating of 
the substance and acceleration of the respective reaction)^. 

The most interesting example of the possible existence of such fluctuations 
is the trace of power like behaviour of the transverse momentum distribution 
in multiparticle production processes encountered in heavy ion collisions 1131a. 
Such collisions are of special interest because they are the only place where 
new state of matter, the Quark Gluon Plasma, can be produced lla. Trans- 
verse momentum distributions are believed to provide information on the 
temperature T of reaction, which is given by the inverse slope of dN/dpT, 
if it is exponential one. If it is not, question arises what we are really mea- 
suring. One explanation is the possible flow of the matter, the other, which 
we shall follow here, is the nonextensivity.-(t>r rather fluctuations leading to 
it). Namely, as was discussed in detail in 113 the extreme conditions of high 
density and temperature occuring in ultrarelativistic heavy ion collisions can 
invalidate the usual BG approach and lead to q > I, i.e., to 



dN(p T ) 

— = const • 

dpT 



T 



(15) 



e It should be noticed that in the case of q < 1 temperature does not reach stationary state 
because, cf. Eq. (0), (1/T) = 1/(T - e/ot'), whereas for q > 1 we had < 1/T >= 1/T . 
As a consequence the corresponding Levy distributions are defined only for e G (0,Too/) 
because for e — * Too', < T >^ 0. Such asymptotic (i.e., for t/r — > 00) cooling of the 
system (T — > 0) can be also deduced form Eq. (h4f) for e —* T^ct' . 
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T = const 



(a) 



a T = 250 MeV 
O T = 200 MeV 

o!l5 V0 1.5 

P T [GeV] 



T = 200 MeV 



(b) 



q = 1.05 
q = 1.0 

V0 1.5 

P T [GeV] 



Figure 1. (a) Normal exponential pj< distributions (i.e., q = 1) for T = 200 MeV (black 
symbols) and T = 250 MeV open symbols), (b) Typical event from central Pb + Pb at 
Ebeam = 3 A-TeV (cf. text for other details) for = 200 MeV for q = 1 (black symbols) 
exponential dependence and q = 1.05 (open symbols). 



Here m is the mass of produced particle and T is, for the q = 1 case, the 
temperature of the hadronic system produced. Although very small (\q — 
1| ~ 0.015) this deviation, if interpreted according to eq. (||)), leads to quite 
large relative fluctuations of temperature existing in the nuclear collisions, 
AT/T ~ 0.12. It is important to stress that these are fluctuations existing 
in small parts of hadronic system in respect to the whole system rather than 
f the event-by-event type for which, AT/T = 0.06/^/N — > for large N (cf. 
for relevant references). Such fluctuations are potentially very interesting 
because they provide direct measure of the total heat capacity C of the system, 



(16) 



(/3 = i) in terms of u> = q— 1. Therefore, measuring both the temperature of 
reaction T and (via nonextensivity q ^ 1) its total heat capacity C, one can 
not only check whether an approximate thermodynamics state is formed in a 
single collision but also what are its theromdynamical properties (especially 
in what concerns the existence and type of the possible phase transitions 13). 
_ To observe such fluctuations an event-by-event analysis of data is needed 
Two scenarios must be checked: (a) T is constant in each event but be- 
cause of different initial conditions it fluctuates from event to event and (b) 
T fluctuates in each event around some mean value Tq. Fig. 1 shows typical 



crete 'q: submitted to World Scientific on February 1, 2008 6 



event obtained in simulations performed for central Pb + Pb collisions taking 
place for beam energy equal E\, eam = 3 ^4-TeV in which density of particles 
in central region (defined by rapidity window —1.5 < y < 1.5) is equal to 

= 6000 (this is the usual value given by commonly used event genera- 
tors 0). In case (a) in each event one expects exponential dependence with 
T = T event and possible departure from it would occur only after averaging 
over all events. It would reflect fluctuations originating from different initial 
conditions for each particular collision. This situation is illustrated in Fig. la 
where p T distributions for T = 200 MeV (black symbols) and T = 250 MeV 
(open symbols) are presented. Such values of T correspond to typical uncer- 
tainties in T expected at LHC accelerator at CERN. Notice that both curves 
presented here are straight lines. In case (b) one should observe departure 
from the exponential behaviour already on the single event level and it should 
be fully given by q > 1. It reflects situation when, due to some intrinsically 
dynamical reasons, different parts of a given event can have different tempera- 
tures au. In Fig. lb black symbols represent exponential dependence obtained 
for T — 200 MeV (the same as in Fig. la), open symbols show the power- like 
dependence as given by (15) with the same T and with q = 1.05 (notice that 
the corresponding curve bends slightly upward here). In this typical event we 
have ~ 18000 secondaries, i.e., practically the maximal possible number. No- 
tice that points with highest px correspond already to single particles. As one 
can see, experimental differentiation between these two scenarios will be very 
difficult, although not totally impossible. On the other hand, if successful it 
would be very rewarding - as we have stressed before. 

One should mention at this point that to the same cathegory of fluctuating 
temperature belongs also attempt Eil to fit energy spectra in both the longitu- 
dinal and transverse momenta of particles produced in the e + e~ annihilation 
processes at high energies, novel nonextepsive formulation of Hagedorn sta- 
tistical wpdel of hadronization process Ejllil and description of single particle 
spectra HL3. 



2.3 Nonexponential decays 

Another hint for intrinsic fluctuations operating in the-physical system could 
be the known phenomenon of nonexponential decays Spontaneous decays 
of quantum-mechanical unstable systems cannot be described by the pure 
exponential law (neither for short nor for long times) and survival time prob- 
ability is P{t) oc t~ s instead of exponential one. It turns out LJ that by using 
random matrix approach, such decays can emerge in a natural way from the 
possible fluctuations of parameter 7 = 1/r in the exponential distribution 
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P(t) = cxp(— jt). Namely, in the case of multichannel decays (with v channels 
of equal widths involved) one gets fluctuating widths distributed according to 
gamma function 



i 



r(|) V2 < 7 >J V 2 < 7 > 



exp 



V) 



2 < 7 > 



(17) 



and strength of their fluctuations is given by relative variance ^ 7 < ^"Ja > ' > ^ 



^, which decreases with increasing v. According to @, it means therefore that, 
indeed, 



Lq(t,T ) = 



2-q 

TO 



1 - (1-q)- 



(18) 



with the nonextensivity parameter equal to q = 1 + -. 



3 Summary 

There is steadily growing evidence that some peculiar features observed in 
particle and nuclear physics (including cosmic rays) can be most consistently 
explained in terms of the suitable applications of nonextensive statistic of 
Tsallis. Here^jSEfi were able to show only some selected examples, more can 
be found in ElO. However, there is also some resistance towards this idea, 
the best example of which is provided in E3. It is shown there that mean 
multiplicity of neutral mesons produced in p — p collisions as a function of 
their mass (in the range from m, ( = 0.55 GeV to My = 9.5 GeV) and the 
transverse mass tot spectra of pions (in the range of tot — 1 -r 15 GeV), both 
show a remarkable universal behaviour following over 10 orders of magnitude 
the same power law function C ■ x~ p (with x = to or x = mx) with P ~ 10.1 
and P ~ 9.6, respectively. In this work such a form was just postulated 
whereas it emerges-.naturally in ^-statistics with q = 1 + 1/-P ~ 1-1 (quite 
close to results of c3). We regard it as new, strong imprint of nonextensivity 
present in multiparticle production processes (the origin of which remains, 
however, to be yet discovered). This interpretation is additionally supported 
by the fact that in both cases considered in c3 the constant c is the same. 
Apparently there is no such phenomenon in AA collisions which has simple 
and interesting explanation: in nuclear collisions volume of interaction is much 
bigger what makes the heat capacity C also bigger. This in turn, cf. cq.(|l6|), 
makes q smaller. On should then, indeed, expect that qhadronic >> Qnuciean 
as observed. 
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As closing remark let us point out the alternate way of getting nonexterLj 
sive (i.e., with q ^ 1) distributions for—thermal models (cf. our remarks in u 
and the more recent ideas presented in Ej) . Notice that if we allow for the tem- 
perature T to be energy E dependent, i.e., that T = T(E) =T + a- (E — E) 
(with a = 1/CV) then the usual equation on the probability P(E) that a 
system A (interacting with the heat bath A' with temperature T) has energy 
E, 

d\n{P(E)]~~dE/T P (E)~exp(-lp\ (19) 

becomes 

E-E 



dlnP(E) ^clE P(E) 

T + a{E-E) 



!-(!-<?)- 



To 



(20) 

with q = 1+a. This approach could thcn_find its possible application to studies 
nf fluctuations on event-by-event basis Ej (with all reservations expressed in 
t3 accounted for). 
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